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ABSTRACT 

We present an upper limit on the polarization of the Cosmic Microwave Background at 7° angular 
scales in the frequency band between 26 and 36 GHz, produced by the Polar experiment. The campaign 
produced a map of linear polarization over the R.A. range 112° — 275° at declination 43°. The model- 
independent upper limit on the E-mode polarization component of the CMB at angular scales £ — 2 — 20 
is 10/iK (95% confidence). The corresponding limit for the B-mode is also 10/.tK. Constraining the 
B-mode power to be zero, the 95% confidence limit on E-mode power alone is 8/iK. 

Subject headings: cosmic microwave background: - cosmology: observations - polarization 



1. INTRODUCTION 

The Cosmic Microwave Background (CMB) is completely 
specified by three characteristics: its spectrum, the spatial 
distribution of its total intensity, and the spatial distribu- 
tion of its polarization. In the past decade, measurements 
of its spectrum and temperature anisotropy have ushered 
in an era of precise determinations of numerous cosmo- 
logical parameters. CMB polarization has the potential 
to provide powerful cross-checks as well as to improve the 
accuracy with which cosmological parameters can be mea- 
sured, most strikingly those related to gravitational waves 
and the ionization history of the Universe. Both of these 
phenomena imprint a signature on the polarization of the 
CMB at the largest angular scales £ < 30. We have built a 
polarimeter called Polarization Observations of Large An- 
gular Regions (Polar) which is optimized to study CMB 
polarization at these scales. 

The polarization of the CMB is a unique probe of ra- 
diative transport in the pre-galactic plasma. Thomson 
scattering of anisotropic CMB radiation by free electrons 
produces CMB polarization, which therefore carries in- 
formation about the Universe during the periods when 
it was ionized. Information about the early (z < 1000) 
ionized epoch is encoded in degree-scale polarized fluctua- 
tions, whereas information about reionization is imprinted 
on scales of tens of degrees, corresponding to the horizon 
scale at the time. Therefore, CMB polarization provides 
the crucial evolutionary link from the recombination epoch 
(z ~ 1000) to the period of large scale structure forma- 
tion (z ~ 6). In particular, reionization should produce 
a new polarized peak near £ < 20, the peak location de- 
pending on the redshift Zi at which the first luminous ob- 
jects heated and reionized the Universe (Zaldarriaga 1998; 
Keating et al. 1998). Despite its fundamental importance, 
Zi is still quite poorly constrained observationally: the lack 
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of an observed Lyman-a trough (due to neutral hydrogen) 
in the spectra of distant quasars (Gunn & Peterson 1965) 
places a lower limit z% > 6, whereas the lack of observed 
suppression of the unpolarized acoustic peaks gives an up- 
per limit Zi < 20 (Wang, Tegmark, & Zaldarriaga 2001). 
The purpose of this Letter is to present the strongest limits 
to date on CMB polarization on the angular scales £ < 30 
relevant to these polarized peaks. 

2. instrument 

Polar's design builds on techniques developed in pre- 
vious searches for CMB polarization (Nanos 1979; Lubin 
& Smoot 1980; Wollack et al. 1993) and is driven by the 
magnitude and angular scale of the anticipated CMB sig- 
nal, and rejection of potential systematic effects. The 
polarimeter is a superheterodyne correlation radiometer 
which detects two orthogonal linear polarization states 
in three radiofrequency (RF) bands in the -Ka-band be- 
tween 26-36 GHz. The two polarization states i € {x,y}, 
Ef F (t, v, fa) — Ei cos(2-Kvt + 4>i) enter a single-mode cor- 
rugated feedhorn and are separated by an orthomode trans- 
ducer (OMT). The OMT's high polarization isolation (> 
30 dB) ensures low cross-polarization. Ef- F are amplified 
by separate HEMT amplifiers (Pospieszalski 1992) cooled 
to 25 K by a mechanical cryocooler. Downconversion from 
the RF band to an intermediate frequency (IF) band (2- 
12 GHz) is performed by Schottky diode mixers, driven 
by a Gunn diode local oscillator (LO) at 38 GHz. In 
the IF band the two polarization states are amplified pro- 
ducing E IF oc Ef- ¥ ', and filtered into three separate IF 
bands, denoted Jl, J2, and J3. The IF bands translate 
into RF bands: Jl (32-36 GHz), J2 (29-32 GHz), and 
J3 (26-29 GHz). Prior to filtering, two diode detectors 
measure the total power of each polarization state, which 
monitors atmospheric opacity. E IF and E IF are corre- 
lated by three Schottky Diode analog multiplier circuits. 
The phase of the LO is switched between and n at 1 
KHz prior to mixing the E^ F waveform. The voltage pro- 
duced by the correlators at this stage thus switches be- 
tween kE^ f ' Ey F and —nE^- F Ey F at 1 KHz, where K is 
the intensity-to-voltage conversion factor. Phase-sensitive 
detection of this modulated signal reduces the effects of 
low frequency drifts in the HEMTs and subsequent com- 
ponents to negligible levels. After low-pass filtering, we 



record an audio-band signal with a DC component Idc = 
K(E^ F Ey F ), where the brackets denote a time-average, 
and AC components proportional to the thermal noise 
from the radiometer, atmosphere, and celestial signals. 
These signals are referred to as the "science channels". 
A second lock-in amplifier for each correlator is referenced 
to the same 1 KHz waveform but delayed in phase by ir/2 
with respect to the phase switch. These signals, here- 
after referred to as "quadrature phase channels" (QPC), 
contain only the noise terms of the RF band and no op- 
tical or celestial signals. The QPC are powerful probes 
of systematic effects produced solely by the radiometer 
and post-detection stages. The output from the QPC is 
proportional to the noise equivalent temperature (NET) 
of the ins trumen t. For a correlation radiometer this is 
NET = v /2/Aiy(T Rx + T A „t)[K - s 1 / 2 ], where Av is the 
bandwidth of the radiometer, Trx is the receiver noise tem- 
perature, and TAnt is the antenna temperature of observed 
optical sources, including diffuse sources such as the atmo- 
sphere and the CMB itself. A schematic outline of Polar 
is presented in Keating et al. (1998), and Keating (2000) 
presents the complete design. The specifications of Polar 
are summarized in Table 1. 

Absolute calibration of the polarimeter is accomplished 
by placing a 0.003" polypropylene sheet in the near-field 
of the antenna, tilted at an angle of 45° relative to the 
optical axis. The dielectric sheet introduces a polarized 
signal into the antenna due to the difference in reflection 
between the parallel and perpendicular electric fields. We 
calibrate the polarimeter by rotating it around the optical 
axis. Near-field, absolute calibration of previous CMB po- 
larimctcrs have used wire-grids to reflect the fields parallel 
to the wire axes (Lubin & Smoot 1980). The dielectric 
sheet material and thickness are chosen to produce an ab- 
solute polarized antenna temperature of ~ 500 mK versus 
200 K for wire-grid calibrators. This provides a calibra- 



tion reference near the NET of the radiometer (accurate 
to 10%), which reduces the dynamic range requirements 
on the IF and post-detection stages. A \ 2 minimization 
fit to a radiometer model is computed during calibration, 
which determines not only the radiometer flux scale, but 
also inter-channel cross-talk and gain imbalance between 
the two polarization states (Keating 2000). 

3. OBSERVATIONS 

The Polar campaign collected a total of 746 hours of 
data between 2000 March 11 and 2000 May 29. Both 
Stokes parameters Q and U were observed through the 
same atmospheric column by a single corrugated scalar 
feedhorn cooled to ~ 40K. The feed produces a nearly 
frequency independent beam with a FWHM ~ 7.0°. The 
beam observes the zenith from Madison, Wisconsin, which 
corresponds to a declination of S = 43°01'48". Contami- 
nation by terrestrial radiation is minimized by our choice 
of vertical drift scans and is further mitigated by two lev- 
els of ground shields (one co-rotating with the polarimeter, 
and the other fixed relative to the ground). Thus, radia- 
tion from the ground must diffract over two ground screens 
before reaching the antenna. 

The Stokes parameters completely describe the linear 
polarization state of radiation produced by Thomson scat- 
tering (Chandrasekhar 1960; Hu & White 1997) and are 
modulated by rotation of the entire instrument around its 
optical axis. The radiometer rotates at / = 0.033Hz. The 
Stokes parameters are modulated at 2/ whereas most po- 
tential systematic signals are modulated at / or DC, and 
are easily discriminated from Q and U. Measurement of 
both Q and U is essential for model independent recovery 
of the 'Electric' (or 'Gradient') and 'Magnetic' (or 'Curl') 
components of the polarization power spectrum (Zaldar- 
riaga & Seljak 1997; Kamionkowski, Kosowsky, & Stebbins 
1997). 



Table 1 





Polar Instrument Specifications 


Channel 3 


[GHz] 


Au c 

[GHz] 


FWHM d 
[°] 


T> c 
J pol 

[/*K] 


Q f 
^sky 

[mK - s 1 / 2 ] 


TP-E/H 
J3 
.J 2 
Jl 


31.9 
27.5 
30.5 
34.0 


7.8/8.0 
0.8 
2.0 
1.9 


7.0 

7.7 
7.1 
7.1 


84(28) 
72(14) 
33(11) 


14.0/20.0 
1.9 
1.0 
1.0 



a TP-E and TP-H measure the total power in the E and H polarization 
planes of the Horn/OMT assembly prior to correlation. 

b Channel Band Ccntroid 

c Channcl Bandwidth. 

d E and H-planc Bcamwidths are equal to within 5%. Measured 
fccd/OMT cross-polarization is < — 40dB for all channels. 



=Mean Polarized Offset T, 



pol 



Q + U , where Q and U are the 
Stokes parameter offsets. Numbers in parentheses denote the correspond- 
ing values for the QPC. 

Measured channel NET for a typical clear day with X a -band zenith 
sky temperature T\ tm ~ 15 K. NET measured at Stokes modulation 
frequency 0.065 Hz. 



4. DATA REDUCTION AND MAPS 

The data were sampled at 20 Hz and converted from 
digital units to Stokes Q and U for each 7.5 minute file. 
Each file contains ~ 14 rotations of the polarimeter during 
which the beam moves ~ 1.9°. For each file, the data 
are sorted into discrete angular bins and synchronously 
demodulated with respect to trigonometric functions of 
the rotation angle. In our coordinate system, defined with 
respect to the local geographic coordinates, we perform 
a x 2 minimization of the angular binned data fit to the 
function 

I(0 t ) = I + C cos 9 t + S sin 9 t + Q cos 29 t + U sin 29 t , (1) 

where 9t — 2tt ft. In addition to the Stokes parameters Q 
and U, the terms C and S (which are synchronous with 
the rotation at frequency /) are monitored to determine 
our sensitivity to rotation-synchronous systematic effects, 
and to monitor atmospheric fluctuations. 

Three levels of diagnostics were used to detect and re- 
move contaminated data: 1) housekeeping and weather 
(dewpoint, cloud cover) cuts; 2) time ordered data (TOD) 
cuts; 3) rotation ordered data (ROD) cuts after construct- 
ing the Stokes parameters. Data were collected while the 
sun (moon) was less than 30° (60°) above the horizon. Ap- 
proximately 2% of the total data collected were removed 
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Fig. 1. — Maps of the Stokes parameters vs. R.A. for each fre- 
quency channel in thermodynamic temperature. For display pur- 
poses, the maps are binned into 4° R.A. pixels with Galactic lati- 
tude |6| > 25°, though the analysis uses 2° pixels. The sky coverage 
was determined by our scan strategy and data selection criteria. 
The Stokes parameters are presented in accordance with the IAU 
definition. 

due to work on the receiver or closure of the observatory 
dome. 1.5% of the data were removed due to completely 
overcast cloud cover. Although atmospheric emission is 
not expected to have significant linear polarization at 30 
GHz (Keating et al. 1998), the correlators are sensitive to 
changes in the total antenna temperature. This coupling is 
caused by imperfect isolation between the E and H planes 
of the OMT and results in a weak dependence of the corre- 
lators on their total bias power (Carretti et al. 2001). The 
bias power is proportional to the system temperature, in- 
cluding atmospheric emission. Since the atmospheric fluc- 
tuation spectrum at 30 GHz falls approximately as 1//" 
(a > 1), due to fluctuations in the amount of prccipitablc 
water vapor, there is residual power at the Stokes param- 
eter modulation frequency 0.065 Hz. This effect is corre- 
lated between all frequency channels, and produces a slow 
drift in the demodulated Stokes parameters, as well as in 
the rotation synchronous terms S and C in equation 1. 
We require that the rotation synchronous terms S and C 
be less than 2xNET. This cut (applied after the TOD has 
been assembled into the ROD) removes almost all weather- 
related contamination of the science channels (35% of the 
total data). Since the basis functions for S and C are or- 
thogonal to those of Q and U, we are assured that only 
data collected during poor weather, and not polarized ce- 
lestial signals, are flagged. Further data files (2.6%) are 
removed if any samples in the file deviate more than 5a 
from the TOD-mean, or 3cr from ROD-mean. After appli- 
cation of all cuts, 121 hours (16%) of the total collected 
data survive. The science channels display mean inter- 
channel cross-correlation coefficients of 1% in the TOD, 
primarily caused by correlated atmospheric fluctuations or 
HEMT gain fluctuations (Wollack & Pospieszalski 1998). 
The QPC display similar correlations. Correlations be- 
tween Stokes parameters in the ROD average 10%, for the 
science channels and ~ 3% for the QPC. For each section of 
data that passes the data quality tests, the time-ordered 
data are converted to a map and covariance matrix via 
the minimum variance map-making procedure outlined in 
Tegmark (1997). The 25 total sections range in length 
from 4-24 hours. An overall offset is removed from each 
section, and the data covariance matrix is adjusted ac- 



cordingly (Tegmark 1997; Bond, Jaffe, & Knox 1998). All 
sections are combined into a final map and covariance ma- 
trix for each channel. Figure 1 presents the maps of the 
Stokes parameters for all correlator channels versus RA. 

5. FOREGROUNDS 

The observations described were conducted over a wide 
range of Galactic latitudes and therefore there is a poten- 
tial for galactic contamination, especially at low latitudes 
(Bennett et al. 1996; Davies et al. 1996). Galactic syn- 
chrotron emission can be up to 75% polarized (Rybicki & 
Lightman 1979). No maps of polarized synchrotron emis- 
sion exist at 30 GHz, and extrapolation of measurements 
at lower frequencies, e.g. Brouw & Spoelstra (1976), is 
not a reliable probe of synchrotron polarization at 30 GHz 
due to Faraday rotation. Although the unpolarized inten- 
sity is not necessarily correlated with the polarized inten- 
sity, we limit our susceptibility to polarized synchrotron 
emission by only using data corresponding to Galactic lat- 
itudes \b\ > 25°, i.e. regions of low unpolarized intensity. 
For reference, extrapolation of the 408 MHz Haslam map 
assuming 10% polarization over our observing region with 
b > 25° yields an RMS of only a couple of /iK. This is 
consistent with our null result described below, indicating 
that we see no evidence of foreground contamination. 

6. LIKELIHOOD ANALYSIS FOR POLARIZED CMB 
FLUCTUATIONS 

A Baycsian maximum likelihood analysis with a uniform 
prior is used to test for polarized fluctuations. Since both 
Q and U are measured simultaneously, we obtain joint 
likelihoods for E and B in a model independent fashion 
(Zaldarriaga 1998). The data from the three frequency 
channels are combined into a single map for each Stokes 
parameter using the minimum variance map-making pro- 
cedure which includes all noise correlations. We form 
a 2iVpi x element data vector x^ = (Qi,Ui) where i = 
{1 . . . Npix}. The likelihood of the model given the data 
is L oc exp(— x T C~ 1 x/2)/|C| 1 / 2 . The covariance matrix 
Cy = Sy + Ny, where Ny = (x^Xj) is the noise covari- 
ance matrix which encodes correlations for both Stokes pa- 
rameters, at each pixel. The noise covariance matrix also 
accounts for the offset removal required to combine maps 
made over long timescales. Ny is a 2A p j X x 2Apj x symmet- 
ric matrix, with 7V p ; x = 84 pixels of width 2° for the |6| > 
25° Galactic latitude-cut data used for the CMB analy- 
sis. Sy is the signal covariance matrix constructed from 
a simple flat band power model parametrized by Te,Tb- 
To place limits on the polarization we consider spectra 
for both E and B of the form £(t + \)Cf /2ir = T\ with 
X G {E, B}. The Q and U signal correlations are com- 
puted by transforming the data coordinate basis to a local 
basis for each pair of pixels. The local basis exploits the 
intrinsic symmetries of the Stokes parameter correlation 
functions (Zaldarriaga 1998; Kamionkowski, Kosowsky, & 
Stebbins 1997). For each pixel pair the Stokes parameters 
in the global basis are rotated such that the axis defining 
the positive Q direction of each pixel lies along the great- 
circle connecting the pixels. The theory covariance matrix 
for all pairs of Stokes parameters depends only on the sky 
coverage and the underlying E and B-mode power spectra: 



(xiXj)-R( ai j)M(r i .? j )R(a ji ) a 



(2) 



Here M is the covariance matrix using a (Q, C/)-convcntion 

where the reference direction is the great circle connecting 

the two points, and R(a) are the rotation matrices defined 

/ cos 2a sin 2a \ 
by R(a) = I . o which rotate the Stokes 

parameters of the i th pixel into a global reference frame 
where the reference directions are meridians (Tegmark & 
de Oliveira-Costa 2000). Sy is constructed from the indi- 
vidual 2x2 blocks of equation (2) by looping over all pixel 
pairs, which requires -^ D i x (-^nix + -0/2 operations. The 
M-matrix is given by (Zaldarriaga 1998) 
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z = r i ■ Yj is the cosine of the angle between the two pixels 
under consideration, Bg = exp[— £{£ + l)a B /2], ob is the 
beam dispersion = 0.425 x FWHM, and Fi,F 2 are func- 
tions of Lcgcndrc polynomials as defined in Zaldarriaga 
(1998) and Tegmark & de Oliveira-Costa (2000). For each 
point in the Tg — Tb plane, equations 4 and 5 are eval- 
uated and substituted into equations 3 and 2 to compute 
Sy. The joint likelihood for {Te,Tb} is shown in Figure 
2. The 95% confidence limits are obtained by integrating 
the likelihood function to determine the region which con- 
tains 95% of the probability. The limits are Te < 10/iK 
and Tb < 10/iK at 95% confidence for the multipole range 
2 < £ < 20. These limits do not include the (~ 10%) cali- 
bration uncertainty. Also shown in figure 2 are the corre- 
sponding limits obtained from our null channels; the QPC. 
Finally, since the B-mode power is expected to be sub- 
dominant, even at large angular scales with no reioniza- 
tion, we have also calculated the limits on Te with Tb = 0. 
Integrating the 1-D likelihood curve for Te with Tb = 
we find Te < 8/iK at 95% confidence. 

7. DISCUSSION 

Polar probes the power spectra at large angular scales 
2 < £ < 20 where the signature of reionization and gravita- 
tional waves are predicted to be most pronounced. While 
we find no evidence for detection of either phenomenon, 
Polar' S limits are the most restrictive upper limits at 
these angular scales, and represent the first limits on E 
and B-modes from an experiment simultaneously mea- 
suring both Q and U. In January 2001 Polar was up- 
graded to measure CMB polarization at 20' scales. The 
upgrade, called COsmic Microwave Polarization at Small 
Scales (Compass), uses the same polarimeter as Polar 
and is an exciting complement to the results presented 
here. Assuming similar per-pixel sensitivity to Polar and 
a generic E-mode polarization spectrum, Compass has the 
potential to detect polarization of the CMB and dramati- 
cally enhance the scientific returns of this nascent field. 
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Fig. 2. — Normalized likelihood contour plots in the Te — Tb 
plane. The countours enclose are 68%, 95.4%, and 99.7% of the 
total probability, corresponding to 1,2, and 3 standard deviation 
intervals, as labeled. The solid lines correspond to the joint like- 
lihood for the combined three frequency channels, and the dashed 
lines are the corresponding null-channel (QPC) likelihood. 
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